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Abstract 

Starting with the relative entropy for two close statistical states we 
define the metric and calculate the scalar curvature R for systems with 
classical, boson and fermion fractal distribution functions with moment 
order parameter q. In particular, we find that for q ^ 1 the scalar cur- 
vature is closer to zero implying that the fractal bosonic and fermionic 
systems are more stable than the standard ones. 
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1 Introduction 

From the theory of fractals pQ we learned that given a statistical weight ft(q, 5) 
of a system with order parameter q and resolution 8, the fractal dimension is 
defined as the exponent d = D q which will make the product lim$->o£l(q, 5)5 d 
finite. With use of the definition of the Boltzmann entropy S(q,5) — lnf2(q, S), 
the relation between the entropy and the fractal dimension D q is given by 

D q = - hm — j — — . (1) 
d^o In 6 

Based on these definitions and with use of the Boltzmann's H theorem, the 
generalized entropy and distribution functions for classical and quantum gases 
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were calculated in Ref. [3]. The average number of particles with energy e was 
shown to be given by 



where a — for the classical case, and the values a = — 1 and a — 1 correspond 
to Bose-Einstein and Fermi-Dirac cases, respectively. For q = 1, Equation ([2]) 
becomes the standard textbook result for classical and quantum ideal gases. The 
distribution functions in Equation ([2J were also obtained in Ref. [3] by consider- 
ing a dilute gas approximation to the partition function of a non-extensive sta- 
tistical mechanics originally proposed in Ref. [3] . It is our purpose to study some 
of the geometric properties of systems with average particle number according 
to Equation © . The idea of using geometry in thermodynamics is not new [5]- 
[9], and several authors developed formalisms to measure the distance between 
equilibrium states through the definition of a metric and the calculation of the 
corresponding scalar curvature as a measure of the interactions [10]- [20]. Some 
of the applications include classical and quantum gases [H] p.6 21 [22 , magnetic 
systems [53]- non-extensive statistical thermodynamics [27 28 [29 , anyon 
gas [5D] [S], fractional statistics [35] and deformed boson and fermion systems 
[33] . Some of the basic results of these approaches include the relationships 
between the metric with the correlations of the stochastic variables, and the 
scalar curvature R with the stability of the system, and the facts that the scalar 
curvature R vanishes for the classical ideal gas, R > 0(R < 0) for a boson 
(fermion) ideal gas, and it is singular at a critical point. Here, we wish to study 
systems with an average particle number given in Equation ([2]) . In Section [2] 
we briefly describe the formalism of systems with fractal distribution functions 
as reported in Ref. [34] . In Section [3J we obtain the metric from the second 
order term in the expansion of the relative entropy between two close statistical 
states, and in Section 3] we use the metric to compute the scalar curvature for 
the classical, Bose-Einstein and Fermi-Dirac cases for q ^ 1. In Section [5] we 
summarize our results. 

2 Fractal models 

2.1 Classical case 

We will use the short notation 



< n(e) > 



1 



(2) 



[l + y %-l)(e-/i)]V(<r-i)+ a ' 



p l = [l + (q-l)P(e l -n)] 1 ^- 1 l 



(3) 



The probability density is defined 



1 




(4) 



P 



Zmb 
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where the partition function 

z MB = n e 



1=0 ni=0 



n« 

1=0 



From the definition of the average number of particles with energy t\ 

Hi ! 

ni Pl 



< ni >-- 



we find after summing the series that 



< ni >= p; 



(5) 



(6) 



(7) 



as required. In the thermodynamic limit we write for the average total number 
of particles < N > 

3/2 ,,oo 



< N > -- 



4ttV ( 2m 



h 3 \P(q-l) 
leading to the expression 



x dx 



[l + X 2 -^-!)^] 1 /^)' 



< N >- 



-2nV 



2m 



3/2 



P(q-l)J [1 - (q - l)^/i] 1 /(9-l)-(3/2) 



5'. 



where S is a series independent of /3^, given by 



S = C + Y, 



(-iy 



in 



m 

2 



C I, 



(8) 



0) 



(10) 



with C m = 



-i/( g _i)_ m -i-(3/2) ■ Solving Equation (|8]) we find that the fugacity 
z — has a temperature dependence given by 



In z = 



1 



1 



1 - 



-2irV 



2m 



h 3 U?-l)/3 



3/2 



s 



(11) 



where oj — 2 (q-i) • F rom Equation (fTTj) we see that the fugacity is restricted to 
the interval < z < e 1 ^ q ~ 1 \ which serves as a cut-off that avoids a negative 
average occupation number. The correct definition of the average energy is 
given by 

4ttV 



< e>- 



00 
2m 



< n(jp) > q p 2 dp, 



(12) 



leading, with Equation ©, to the required classical result <e>=|<iV>fcT 
[33] . We should remark, that Equations (O and (fT2")l are related by the standard 
definition 

9 In Z M b 

<£>= dP~- (13) 
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2.2 Boson case 

Similarly to the classical case, the definition 



— n. 

^ V n 3 

j=0 Z^n j= P3 



(14) 



leads to the average occupation number 



< n j >= ~r—r (15) 
Pi - 1 



with the probability density and the partition function 



P 



Zbe . , 

3=0 



j=0 F 3 

As in the standard, 5 = 1, Bose-Einstein case the chemical potential is negative. 



2.3 Fermion case 

For the Fermi-Dirac case the average occupation number 

1 



pi + 1 



(18) 



is obtained by defining 



P = 



-\\^: ■ as) 



Zfd . ,, 

3=0 

z *° = liter 

j=0 n j= Q 

= n( i+ ^ Ti )' ( 2 °) 

3=0 

and the requirement that the average occupation number < nj ><G [0, 1] leads 
to restrict the fugacity to the interval < z < e l ^ q ~ x \ For Bose-Einstein and 
Fermi-Dirac cases we obtain 

-^=->+i±(V H> s+ * < 2i > 

' 1=0 k=l v 7 



1 We should remark that our partition function differs from that reported in Ref . [3] . 
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3 The metric 



For the three cases discussed in the previous section we cannot adopt any of the 
standard definitions for the metric like for example |21| 

<9 2 lnZ , , 
^^Wd^ 1 Z? 1 = /? 5 /3 2 = -/3M (22) 

which is valid for exponential distributions. 

The relative entropy, H(p\\P) = J2iPi (~^ n (^")) > * s a ver y useful concept. 
For example, for two close distribution functions p(x) and p{x + A) we can 
obtain a Fisher's information measure 

as part of the second order term in A for the entropic form, with p, = 1, 
S = —J2iPi^ n Pi an( i & = l n Pi) A ' ES] where \x is a fractional 

parameter. By defining 0(cc) = p^f and considering the Fisher information 
as a lagrangian density leads to linear and nonlinear differential equations for 
p, = 1 and /i 7^ 1, respectively. 

Here, based on work in Ref.|10) we expand the relative entropy for p = 1 
between two close densities p(f3) and p(j3 + df3) up to second order in df5 a . 
Therefore, the information distance I(p(/3), p[fi + d/3)) between the two close 
states is written 

I(p{n,P{P a + d(3 a )) = Trp (\np(f3 a ) - \np((3 a + d(3 a )) , (24) 
such that expanding the second order term gives for the metric 



WWi+ Trp^<n l > w ^. (25) 



1=0 

A simple inspection shows that in the q — ¥ 1 limit Equation (|25|) reduces to 
Equation (j2"2"j) . Equation (|2"5"j) can be simplified leading to the three correspond- 
ing metrics: 

MB _ \ ^ 1 d Pi d Pi 



(26) 



1=0 

dpi dpi 

P i{pi-iydp a 9/37' 



(29) 



which can be summarized in the general formula 

<ni> 2 dpi dp, 
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Writing in general Pl = [l + (q- 1) J2 a /J ^ ] 1 /^ 1 ' we find 
d d 



d/3 a d/3 x 



InZ = pV^KFi (q<n l >~a< m > 2 ) (31) 



1=0 



g a x = Y.Pi~ 2qF i aF i X ( <n i > -a<m > 2 ) (32) 



1=0 



4 Scalar curvature 

4.1 Classical case 

It has been shown that the scalar curvature [H][IB] vanishes for the standard 
case, but it is tempting to speculate whether that is also the case for q ^ 1 . In 
the thermodynamic limit,with x = f3e, we write for example 

2 f°° T 5 / 2 dr 

3ii = 4wr 2 A- 3 / x —™ — (33) 

Jo [I + (q - l)(x + ~f)} — 



'TV 



where hereafter 7 = — /3/i. This integral converges for 2 ( g -i) - > 0' restricting the 
values of q to the interval q € [1,5/3). With use of the integral representation 
of the T-function 

/>oo 

T(y) = w y t y - l e- wt dt, y > ; w > 0, (34) 
Jo 

we obtain for the components of the metric tensor 



2- 



.912 = Vpr^hz,*, (35) 
522 = V\~ 3 h 1/2 , 



where the function 



2 r(A + l)r(^y — A) 1 

hx= Jt(q-l)W F(E) [1 + (g _ 1)7] ^r-A' ^ 

satisfies 

* =-M,_, (37) 
As is well known |17j . the scalar curvature is given by 

R=^-Ri2i2, (38) 

where cieig = 511522 — 312512 and the non- vanishing part of the curvature tensor 
Rap-yX is given in terms of the Christoffel symbols 



6 




Figure 1: The scalar curvature R, in units of A 3 T^ _1 , as a function of the fugacity 
z for bosons at constant f3 for the cases of q = 1 (solid line), q = 1.1 (dashed 
line) and q = 1.2 (dotted line). 



Rap-yX = {FnaX^ep-f ~ ^rjay^epx) ■ (39) 

A simple calculation leads to the result 

R = m 7 \2 \ ^> h l/2 h l/2 ~ Gh l/2 h 5/2 + ^3/2^-1/2^5/2 ) ) (40) 



4 ( detg )2 ( 5/l l/2^3/2 - 6K/2 h 5/2 + ^3/2^-1/2^5/1 

such that after replacement of the dehnition of the function h\ we get that the 
scalar curvature for the classical fractal case is identically equal to zero. There- 
fore, in this case the parameter q does not play any role as far as correlations 
are concerned. 

4.2 Boson and fermion cases 

Here, in order to evaluate the corresponding scalar curvatures we need to replace 
the summations in Equations (|2"T1) and (|2"8|) by integrals 

3-2q 

Gf = 4= / d*. (41) 

V^io (r!^r±i) 2 

where the + sign is for fermions and the — sign for bosons, and the function 
SI = 1 + (q — l)(x + 7). In particular, the metric component gn is written 

g n = V\- 3 (3- 2 G± 2 , (42) 

and its integral converges for 1 < q < 5/3. 



The functions Gf also satisfy 

dG 



± 



dt = - XG ^ (43) 
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Figure 2: The scalar curvature R, in units of A 3 1^ _1 , as a function of the fugacity 
z for fermions at constant /3 for the cases of q = 1 (solid line), q — 1.1 (dashed 
line) and q = 1.2 (dotted line). 

and thus the corresponding equations for R are equivalent to Equation (|40[) with 
the replacement of the function h\ by the functions G^. Figures 1 and 2 show 
the results of a numerical calculation of the scalar curvature R as a function of 
the fugacity z for the parameter values q = 1,1.1,1.2 for boson and fermions 
respectively. 

5 Conclusions 

In this paper, starting from the relative entropy for two close statistical states 
we defined the metric for systems with fractal distribution functions with order 
parameter q. We calculated the scalar curvature R and found that it vanishes 
for the classical ideal gas, as in the standard case. Numerical calculations for the 
boson and fermion systems show that the corresponding values of R as a function 
of the fugacity z are closer to zero than those in the q — 1 case, implying that 
the departure from the value q — 1 makes the systems more stable. Therefore, 
for q 7^ 1 bosons will be less attractive and fermions less repulsive that their 
standard counterparts. Our results are in agreement with those obtained in a 
previous work [37] wherein we showed that long-range correlations for the fractal 
Bose case decrease when the parameter q departs from the standard value q = 1 . 
On the other hand, if one wishes to consider the order parameter q as a non- 
extensive parameter it has to be within the context of considering these fractal 
systems as a dilute approximation to non-extensive statistical mechanics, which 
consists in replacing the Tsallis partition function by a factorized one. This type 
of approximation has been shown [35] to be good outside a temperature interval 
that shifts to higher values of T when the number of energy levels increases. 
Our results also show that the sign of R remains unchanged as a function of z 
implying that these systems do not exhibit anyonic behavior, a fact that looks 
impossible to check by performing an expansion for z ss to obtain the second 
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virial coefficient because the partition function is a function of lnz. In addition, 
our results contrast with the cases of systems with quantum group symmetry 
where the parameter q interpolates between bosons and fermions in two and 
three dimensions |39| . 

Acknowledgements 

I am grateful to the anonymous reviewers for their comments and constructive 
criticism to improve the original manuscript. 

References 

[1] J. Feder, Fractals (Plenum, New York, 1988). 

[2] F. Buyukkilic and D. Demirhan, Phys. Lett. A 181 (1993) 24. 

[3] F. Buyukkilic, D. Demirhan and A. Giilec, Phys. Lett. A 197 (1995) 209. 

[4] C. Tsallis, J. Stat. Phys. 52 (1988) 479. 

[5] L. Tisza, Generalized Themodynamics (MIT, Cambridge, 1966). 

[6] R. B. Griffiths and J. C. Wheeler, Phys. Rev. A 2 (1970) 1047. 

[7] F. Weinhold, J. Chem. Phys. 63 (1975) 2479. 

[8] S.-I. Amari, Differential- Geometrical Methods in Statistics (Springer- 
Verlag, Berlin, 1985). 

[9] S.-I. Amari and H. Nagaoka, Methods of Information Geometry (AMS, 
Rhode Island, 2000). 

[10] R. S. Ingarden, H. Janyszek, A. Kossakowski and T. Kawaguchi, Tensor 
N.S. 37 (1982) 105. 

[11] G. Ruppeiner, Phys. Rev. A 20 (1979) 1608. 

[12] W. K. Wootters, Phys. Rev. D 23 (1981) 357. 

[13] R. Gilmore, Phys. Rev. A 30 (1984) 1994. 

[14] G. Ruppeiner, Phys. Rev. A 32 (1985) 3141. 

[15] R. Gilmore, Phys. Rev. A 32 (1985) 3144. 

[16] J. Nulton, P. Salamon, Phys. Rev. A 31 (1985) 2520. 

[17] S. Weinberg, Gravitation and Cosmology (Wiley and Sons, Inc., New York, 
1972). 

[18] H. Janyszek, Rep. Math. Phys. 24 (1986) 1; Rep. Math. Phys. 24 (1986)11. 



9 



[19] H. Janiszek and R Mrugala, Rep. Math. Phys. 27 (1989) 145. 

[20] G. Ruppeiner, Am. J. Phys. 78 (2010) 1170, and references therein. 

[21] H. Janiszek and R Mrugala, J. Phys. A: Math. Theor. 23 (1990) 467. 

[22] D. Brody and D. Hook , J. Phys. A: Math. Theor.42 (2009) 023001. 

[23] H. Janiszek and R Mrugala, Phys. Rev. A 39 (1989) 6515. 

[24] H. Janyszek, J. Phys. A:Math. 23 (1990) 477. 

[25] D. Brody and N. Rivier, Phys. Rev. E 51 (1995) 1006. 

[26] W. Janke, D. A. Johnston and R. Kenna, Physica A 336 (2004) 181. 



[27] R. Trasarti-Battistoni, cond-mat/0203536 



[28] M. Portesi, A. Plastino and F. Pennini, Physica A 365 (2006) 173. 

[29] A. Ohara, Phys. Lett. A 370 (2007) 184. 

[30] B. Mirza and H. Mohammadzadeh, Phys. Rev. E 79 (2008) 021127. 

[31] B. Mirza and H. Mohammadzadeh, Phys. Rev. E 80 (2009) 011132. 

[32] B. Mirza and H. Mohammadzadeh, Phys. Rev. E 82 (2010) 031137. 

[33] B. Mirza and H. Mohammadzadeh ,J. Phys. A: Math. Theor. 44 (2011) 
475003 

[34] M. R. Ubriaco, Phys. Rev. E 60 (1999) 165. 

[35] G.V Vstovsky, Phys. Rev. E 51 (1995) 975. 

[36] M. R. Ubriaco, Phys. Lett. A 373 (2009) 4017. 

[37] M. R. Ubriaco, Phys. Rev. E 62 (2000) 328. 

[38] Q.-A. Wang and A. Le Mehaute, Phys. Lett. A 235 (1997) 222. 

[39] M. R. Ubriaco, Phys. Rev. E 55 (1997) 291. 



10 



